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! Abstract 
(N ■ 

^ : 

O , We propose a new set of supersymmetric benchmark scenarios, taking into account the 

^ . constraints from LEP, b — > s-y, — 2 and cosmology. We work in the context of the 

^ ' constrained MSSM (CMSSM) with universal soft supersymetry-breaking masses and assume 
that R parity is conserved. We propose benchmark points that exemplify the different 
D \ generic possibilities, including focus-point models, points where coannihilation effects on the 
^ . relic density are important, and points with rapid relic annihilation via direct-channel Higgs 
' poles. We discuss the principal decays and signatures of the different classes of benchmark 
\-{ ' scenarios, and make initial estimates of the physics reaches of different accelerators, including 
' ' ' the Tevatron collider, the LHC, and e~^e colliders in the sub- and multi-TeV ranges. We 
stress the complementarity of hadron and lepton colliders, with the latter favoured for non- 
strongly-interacting particles and precision measurements. We mention features that could 
usefully be included in future versions of supersymmetric event generators. 
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1 Introduction 



The completion of the LEP experimental programme brings to an end an era of precise 
electroweak measurements and the search for new particles with masses ^ 100 GeV. With 
the start of Run II of the Fermilab Tevatron collider, the advent of the LHC and perhaps 
a linear e~^e~ collider, the experimental exploration of the TeV energy scale is beginning in 
earnest. 

The best-motivated scenario for new physics at the TeV energy scale is generally agreed to 
be supersymmetry. Theoretically, it is compellingly elegant, offers the possibility of unifying 
the fermionic matter particles with the bosonic force particles, is the only framework thought 
to be capable of connecting gravity with the other interactions, and appears essential for the 
consistency of string theory. However, none of these fundamental arguments offer clear 
advice as to the energy scale at which supersymmetric particles might appear. 

The first such argument was provided by the hierarchy problem: if supersymmetric par- 
ticles weigh less than of order 1 TeV, they may stabilize the electroweak scale mz <^ mp 
10^^ GeV. The heavier the supersymmetric particles, the more the fine-tuning of the model 
parameters required to fix mz at its observed value. However, it is difficult to attach quan- 
titative significance to any specific measure of the amount of this fine tuning. 

LEP has provided no direct evidence for any physics beyond the Standard Model, but it 
has provided several indirect hints that supersymmetry may indeed lie around the corner. 
One such hint was provided by LEP's very accurate measurements of the gauge couplings, 
which are highly consistent with a supersymmetric Grand Unified theory (GUT) if the su- 
persymmetric partners of the Standard Model particles weigh less than about 1 TeV as 
suggested by the hierarchy problem p[. Secondly, the precise electroweak data from LEP 
and elsewhere suggest that the Higgs boson is relatively light 0: 

mn = 98+1 GeV (1) 

in good agreement with the prediction of the minimal supersymmetric extension of the Stan- 
dard Model (MSSM), if the squarks weigh ^ 1 TeV. Direct searches at LEP provided the 
lower limit ttih > 113.5 GeV (95% CL) In the final weeks of its run, LEP provided 
tentative evidence for a Higgs boson weighing 115.0lo!9 GeV [Q, perfectly consistent with 
the range (|l]) expected on the basis of the precise LEP electroweak measurements, as well as 
with the MSSM. Indeed, the effective potential of the Standard Model would be so sensitive 
to destabilization by radiative corrections if this tentative LEP evidence were to be con- 
firmed, that some form of supersymmetry would probably be needed to stabilize our familiar 
electroweak vacuum ||^. 
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A completely independent motivation for supersymmetry at the TeV scale is provided 
by the cold dark matter advocated by astrophysicists and cosmologists. If R parity is con- 
served, as we assume here, the lightest supersymmetric particle (LSP) is an ideal candidate 
to constitute the cold dark matter, if it weighs ^ 1 TeV. We assume here that the LSP is 
the lightest neutralino x [§!• The relic LSP density increases with the masses of the super- 
symmetric particles, so the cosmological upper limit Q^h"^ < 0.3 may, in principle be used to 
set an upper limit on the sparticle masses. However, in practice, one must be careful not to 
discard 'funnels' in the MSSM parameter space where heavier sparticles may be permitted. 

Finally, we should add that the recent precise measurement of the anomalous magnetic 
moment of the muon, — 2, which is in apparent disagreement with the Standard Model 
at the 2.6-cr level 0, has led to many speculations about new physics at the TeV scale. 
Prominent among these have been various supersymmetric interpretations of the possible 
discrepancy. These offer further encouragement that supersymmetry might be discovered 
at the LHC or before. However, caution advises us to await confirmation of the initial 
experimental value of and to seek consensus on the calculation of hadronic contributions 
to — 2 before jumping to any conclusions. 

Nevertheless, the front-running nature of the supersymmetric candidacy for new physics 
beyond the Standard Model has motivated many studies of its experimental signatures at 
future colliders. In order to focus these discussions, and to provide standards of comparison 
for different analyses, experiments and accelerators, specific benchmark choices of super- 
symmetric parameters have often been proposed. For example, several years ago, several 
such benchmark scenarios were used to evaluate the capabilities of the LHC for detecting 
supersymmetry More recently, analogous benchmarks have been used in linear 

collider studies fill . 

Unfortunately, time has overtaken some of these benchmark scenarios, whose parameters 
have by now been excluded by direct experimental searches for superymmetry and the Higgs 
boson at LEP, or because they predict unacceptable values of (7^ — 2 or 6 ^ 37 decay, or 
because they predict unacceptable values of the relic LSP density VL^h'^. 

The purpose of this paper is to propose a new set of benchmark supersymmetric model 
parameters that are consistent with the experimental constraints, as well as cosmology. They 
may therefore provide helpful aids for understanding better the complementarity of different 
accelerators in the TeV energy range. We restrict our attention to a constrained version 
of the MSSM (CMSSM) which incorporates a minimal supergravity (mSUGRA)-inspired 
model of soft supersymmetry breaking. In the CMSSM, universal gaugino masses mi/2, 
scalar masses mo (including those of the Higgs multiplets) and trilinear supersymmetry 
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breaking parameters Aq are input at the supersymmetric grand unification scale. In this 
framework, the Higgs mixing parameter fi can be derived (up to a sign) from the other 
MSSM parameters by imposing the electroweak vacuum conditions for any given value of 
tan/3. Thus, given the set of input parameters determined by {mi/2, rriQ, Aq, tan P, sgn{fi)}, 
the entire spectrum of sparticles can be derived. Here we will further restrict our attention to 
Aq = 0, for simplicity. We do not consider benchmarks for models with gauge- [|l^], gaugino- 
131 or anomaly- mediated |14] supersymmetry breaking, or for models with broken R parity. 



Studies of these and other models would be interesting complements to this work, and we 
comment on them in the last Section of this paper. 

Fig. |l| illustrates qualitatively the CMSSM benchmark points we propose, superimposed 
on the regions of the (mi/2,mo) plane favoured by LEP limits, particularly on ruh, b ^ sj 
and cosmology. Electroweak symmetry breaking is not possible in the dark-shaded triangular 
region in the top left corner, and the lightest supersymmetric particle would be charged in 
the bottom right dark-shaded triangular region. The experimental constraints on and 
6 ^ 57 exert pressures from the left, as indicated, which depend on the value of tan (3 and 
the sign of fi. The indication of a deviation from the Standard Model in — 2 disfavours 
/i < and large values of mo and mi/2 for /i > 0. The region where Q^h"^ falls within the 
preferred range is indicated in light shading, its exact shape being dependent on the value of 
tan/3, and to some extent on the Standard Model inputs m^, m^, and a^. As discussed later 
in more detail, in addition to the 'bulk' region at low mo and mi/2, there is a coannihilation 
'tail' extending to large mi/2 [TE, TB\, a 'focus-point' region at large mo near the boundary 
of the region with proper electroweak symmetry breaking [|l^ , and narrow rapid-annilation 



'funnels' at intermediate mo/mi/2 for large tan/? [|T8| [T9|, The interplays of 



these features for different values of tan/3, sgn{fi), together with the corresponding proposed 
benchmark points, are shown in Figs. ^ and |^. 

It is possible to learn much from general theoretical scans of the CMSSM parameter 
space, but the LHC experience also showed the complementary advantages of devoting some 
experimental attention to specific benchmark points [§], where the nitty-gritty problems 
of disentangling possible overlapping experimental signals and extracting measurements of 
CMSSM parameters can be confronted. We do not propose here a 'fair' statistical sampling 
of the allowed CMSSM parameter space. Rather, we propose benchmark points that span the 
essential range of theoretical possibilities, given our present knowledge. Some of the points 
we propose might soon become obsolete, for example because of Tevatron data or reductions 
in the error in g^ — 2. As seen in Fig. |I], many of the points we propose are spread over 
the allowed part of the 'bulk' region, at different values of tan jS. However, we also propose 
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mi/2 

Figure 1: Qualitative overview of the locations of our proposed benchmark points in a generic 
(mi/2, ?^o) plane. The light (turquoise) shaded area is the cosmologically preferred region with 
0.1 < ^x^"^ — '-'•3; whose exact shape depends on the value of tan P, and to some extent on 
the Standard Model inputs rrit, mi, and a^. In the dark (brick red) shaded region at bottom 
right, the LSP is the charged fi, so this region is excluded. Electroweak symmetry breaking is 
not possible in the dark (pink) shaded region at top left. The LEP experimental constraints, 
in particular that on rrih, and measurements ofb-^S'y exert pressure from the left side. The 
BNL E821 measurement of — 2 favours relatively low values of rriQ and mi/2 for fi > 0. 
The CMSSM benchmark points we propose are indicated roughly by the (blue) crosses. We 
propose points in the 'bulk' region at bottom left, along the coannihilation 'tail' extending to 
larger mi/2, in the 'focus-point' region at large rriQ, and in the rapid- annilation 'funnel' that 
may appear at intermediate mQ/mi/2 for large tan/3. 
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some points strung along the coannihilation 'tail', including one at the extreme tip, and two 
points each in the 'focus-point' region and the rapid-annihilation 'funnels'. Some of these 



points might appear disfavoured by fine-tuning arguments I^S], but cannot be excluded. 
Taken together, the points we propose serve to highlight the different possible scenarios with 
which future colliders may be confronted. The input parameters for the benchmark points 
we propose, together with the resulting spectra as calculated using the code SSARD ||25|, are 
shown in Table |^. 

The layout of this paper is as follows. In Section 2, we discuss the various experimental 
and other constraints on supersymmetric scenarios, and discuss how we implement them. 
Then, in Section 3, we introduce the set of benchmark scenarios we propose, motivating 
our choices in the multidimensional parameter space of the MSSM. For convenience, we 



introduce in Section 4 versions of these benchmarks calculated with suitable ISASUGRA ^6 
inputs, and we then use ISASUGRA to discuss the decay signatures of heavier sparticles, which 
are quite distinctive in some of these benchmark scenarios. Then, in Section 5 we take first 
looks at the physics reaches of various TeV-scale colliders, including the Tevatron, the LHC, 
a 500-GeV to 1-TeV linear e~^e~ collider such as TESLA, the NLC or the JLC, and a 3- to 
5-TeV linear e~^e~ collider such as CLIC [|. Finally, in Section 6, we review our results on 
the CMSSM benchmark scenarios we propose, discuss some of the future work that might 
be done to investigate further these benchmark supersymmetric scenarios and use them as 
a guide to understanding the physics opportunities offered by future colliders, and mention 
other possible scenarios that could also be studied. 

2 Experimental and Cosmological Constraints 

We implement the experimental and cosmological constraints using a code SSARD that incor- 
porates the two-loop running of the input soft supersymmetry-breaking parameters from the 
input scale Mqut (defined as the scale where gi and g2 meet) down to the electroweak scale, 
identified with rriz- Exact gauge coupling unification is enforced, and the strong couplings 
constant as{mz) is a prediction. The parameter is extracted by minimizing the one-loop 
corrected effective potential p^, ^ at the scale m^, while the pseudoscalar Higgs mass niA 
is computed using the results of ||2^. The radiative corrections to the light Higgs boson 
mass m/j are computed with the FeynHiggs code The full one-loop corrections to the 



physical chargino and neutralino masses are included [pTl p2|, p3[. The code also calculates 



6 ^ S7 |3^, gn — 2 and the cosmological relic density using consistent conventions. We note, 
^We comment in passing on the capabilities of colliders. 
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Table 1: Proposed CMSSM benchmark points and mass spectra (in GeV), as calculated 



using SSARD / |^3| / and FeynHiggs j^/ . The renormalization- group equations are run down 
to the electroweak scale mz, where the one-loop corrected effective potential is computed 
and the CMSSM spectroscopy calculated, including the one loop corrections to the chargino 
and neutralino masses. The pseudoscalar Higgs mass rriA is computed as in J^/ . Exact 
gauge coupling unification is enforced and the prediction for as{mz) is shown (in units of 
0.001). It is also assumed that = and mb{mb)^^ = 4.25 GeV. For most of the points, 
rrit = 175 GeV is used, but for points E and F the lower value mt = 171 GeV is used, for 
better consistency with ^T^J . 
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in particular, that the inclusion of one-loop corrections to chargino and neutralino masses 
is important for implementing accurately the LEP limits and the boundary of the region 
favoured by cosmology. We discuss later the problems encountered in matching the sparticle 
spectra obtained using this and other codes that implement these constraints using different 
approaches and/or approximations. 



2.1 Sparticle Searches 

The most important direct experimental constraints on the MSSM parameter space are 
provided by LEP searches for sparticles |^ and Higgs bosons the latter constraining 
the sparticle spectrum indirectly via radiative corrections, particularly those associated with 
third-generation supermultiplets. We use here the preliminary combined results that are 
based on data-taking at centre-of-mass energies up to about 208 GeV. 



Upper limits at 95% CL on the the cross section for chargino-pair production were set |^ 



for all kinematically accessible chargino masses as a function of the neutralino mass, assuming 
that the branching ratio for W^x^ was 100%. For neutralino masses approximately 

half the chargino mass, the upper limit obtained using 35 pb^^ of integrated luminosity at 
^/s > 207.5GeV is around 0.5 pb. These cross-section limits can be interpreted within the 
MSSM for some specific parameter values; for tan/3 = 2, fi = —200 GeV and sneutrino 
masses exceeding 300 GeV, the lower limit on the chargino mass is 103.5 GeV 0. 

Similarly, the combined LEP data at ^/s from 183 to 208 GeV were used to search for 



sleptons [|3^. Events containing two charged leptons and missing energy were analysed and 
upper limits set on the cross section times branching-ratio squared for slepton-pair production 
followed hj £ i+x decay, as functions of the slepton and neutralino masses. The limits vary 
substantially with the masses, but typically the limits are 40 fb for the selectron and smuon 
search and 100 fb for the stau search. Within the context of the MSSM, these experimental 
limits lead to the exclusion of major portions of the right-handed slepton, neutralino mass 
plane at 95% CL. The mass limits were evaluated for tan/3 = 1.5 and fi = —200 GeV. For a 
neutralino mass of 40 GeV, the lower limits on the right handed slepton masses are 99.4 GeV, 
96.4 GeV and 87.1 GeV for the selectron, smuon and stau respectively 0. 

There are also important constraints on the squark and gluino masses from Run I of the 
Tevatron extending up to 300 GeV if rriq ~ rrig, and additional constraints on stop and 

^ There are also model-dependent limits on the supersymmetric parameter space derived from searches 
for associated xx' production at LEP. 

■^We note that, for both the chargino and the slepton searches, the sensitivity is much reduced for small 
values of the mass difference between the parent sparticle and the neutralino. However, this caveat is 
unimportant for the CMSSM as studied here. 
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sbottom squarks from LEP. Since they do not play a role in our analysis, we do not discuss 
them in detail. 

Analyses indicate that the experimental search results can be interpreted as chargino 
and slepton mass bounds close to the kinematic limits for most of the CMSSM parameter 
range Therefore, we show in panels (a,b) of Fig. |^ the contours (dot-dashed) in the 

(mi/2,'n^o) plane corresponding to m^± = 103.5 GeV and mg = 99 GeV. These contours are 
omitted from the remaining figures, for clarity. 

As mentioned in the Introduction, the lower limit on the mass of a Standard Model 
Higgs boson imposed by the combined LEP experiments is 113.5 GeV @|. This lower limit 
applies also to the MSSM for small tan (3, even if squark mixing is maximal. In the CMSSM, 
maximal mixing is not attained, and the e^e^ —>■ + h production rate is very similar to 
that in the Standard Model, for all values of tan/?. Therefore, the LEP hint for a Higgs 
boson weighing 115.0io!9 GeV 0], which is compatible with the background-only hypothesis 
at the 0.4% CL, can also be interpreted in the CMSSM. 

To calculate theoretically the mass of the lightest MSSM Higgs boson, we use the 
FeynHiggs code [^], which includes one- loop effects and also the leading two-loop con- 
tributions, and gives results that are somewhat higher than those obtained using [4^]. In 



order to account for uncertainties in theoretical calculations of in the MSSM [^| for any 
given value of rrit, we consider this LEP range to be consistent with CMSSM parameter 
choices that yield 113 GeV < rrih < 117 GeV. The theoretical value of rrih in the MSSM 
is quite sensitive to nit, the pole mass of the top quark: we use nit = 175 GeV as default, 
but mention explicitly the cases where rrit = 171 GeV has been used. Calculations of the 
Higgs mass and other quantities are also sensitive to the bottom-quark mass (particularly 
at large tan/5), for which we choose mb{mb)^^^ = 4.25 GeV for the running mass. All but 
one of the benchmark points we propose satisfy rrih > 113 GeV for rrit = 175 GeV. In view 
of the expected accuracy ~ 3 GeV of the FeynHiggs code, we therefore consider that all the 
proposed points are compatible with the LEP lower limit of 113.5 GeV [Q. 

2.2 6 ^ S7 Decay 



We implement ||3^ the new NLO b ^ s'j calculations of when M > 500 GeV, where 



M = Min(mg, rrig). Otherwise, we use only the LO calculations and assign a larger theoretical 
error. For the experimental value, we combine the CLEO measurement with the recent 
BELLE result assuming full correlation between the experimental systematics 0, finding 



^This is conservative, but the available information does not justify a less conservative approach, and this 
assumption is in any case not very important. 
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Figure 2: The (mi/2, mo) planes for tan (5 = (a) 5 (fi > 0), (b) 10 (fi > 0), (c) 10 (fi < 0), 
all for rrit = 175 GeV, and (d) 10 (fi > 0) with rrit = 171 GeV. In each case we have 
assumed Aq = and mb{mb)sM — ^-25 GeV, and used the SSARD code. The near-vertical 
(red) dot-dashed lines are the contours = 113 GeV, as evaluated using the FeynHiggs 
code. The medium (dark green) shaded regions are excluded by b s'j. The light (turquoise) 
shaded areas are the cosmologically preferred regions with 0.1 < n^h^ < 0.3. In the dark 
(brick red) shaded regions, the LSP is the charged fi, so this region is excluded. The regions 
allowed by the E821 measurement of at the 2-a level are shaded (pink) and bounded by 
solid black lines, with dashed lines indicating the 1-a ranges. Electroweak symmetry breaking 
is not possible in the dark (pink) shaded region at the top left of panel (d). The (blue) crosses 
denote the proposed benchmark points A to F. 
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mi/2 (GeV) mi/2 (GeV) 

Figure 3: The (mi/2, mo) planes for tan/3 = (a) 20 (fi > 0), (b) 35 (fi > 0), (c) 35 
(fi < 0), and (d) 50 (fi > 0), found using SSARD and assuming Aq = 0,mt = 175 GeV and 
''^bij^bls^M = 4.25 GeV. The notations are the same as in Fig. ^. The (blue) crosses denote 
the proposed benchmark points G to M. At larger tan /3, the size as well as the exact shape 
of the cosmologically preferred region obtained is subject to considerable uncertainty, and 
different programs yield different answers for the same fixed values of the input parameters. 
The differences arise due to different calculational algorithms, and to neglecting different sets 
of higher-order terms. We elaborate more on these issues in Section |^. 
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B{b S7) = (3.21±0.44±0.26) x 10~^. In our implementation, we allow CMSSM parameter 
choices that, after including the theoretical errors ath due to the scale and model dependences, 
may fall within the 95% confidence level range 2.33 x 10^^ < B{b 57) < 4.15 x 10~^. In 
general, we find in the regions excluded when < that the predicted value of B{b 57) is 
larger than this measured range, whereas, when fi > 0, the exclusion results from B{b —>■ 57) 
being smaller than measured. Table ^ shows the values of B{b S7) calculated in our 
proposed benchmark scenarios. 
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(+ At) 


(292) 


(47) 


(117) 


(177) 


(153) 


(559) 


(97) 


(1098) 


(83) 


(294) 


(537) 


(138) 


(1276) 




6.0 


1.3 


5.7 


7.0 


106 


85 


9.3 


36 


12 


32 


91 


7.3 


33 


(+ At) 


(6.0) 


(1.3) 


(5.9) 


(7.0) 


(372) 


(1089) 


(11) 


(36) 


(13) 


(33) 


(125) 


(29) 


(206) 



Table 2: Derived quantities in the benchmark models proposed. In addition to the relic 
density Q^h"^, the supersymmetric contribution to = [g^ — 2)/2 in units of , and the 
b s'j decay branching ratio 10~^, we also display the amount of electroweak fine-tuning 
( all of the above quantities are calculated using SSARD ), and the amount of electroweak 



fine-tuning, calculated with the BMPZ code j^/ , using the ISASUGRA 7.51 versions of the 
input parameters. 



2.3 Muon Anomalous Magnetic Moment 

The BNL E821 experiment has recently reported |^ a new value for the anomalous magnetic 
moment of the muon: g^—2 = 2xa^, which yields an apparent discrepancy with the Standard 
Model prediction at the level of 2.6 a: 

5a^ = (43 ± 16) X 10"^°. (2) 

The largest contribution to the stated error is due to statistics, and is expected to be re- 
duced soon by a factor two or more. The systematic errors reported by the BNL E821 
experiment are considerably smaller in magnitude. The largest uncertainty in the Standard 
Model prediction is that due to the hadronic contributions: Sa'^"''^ ~ 7 x lO^^''. The largest 
contribution to a'^f''^ is in turn due to vacuum polarization diagrams, with the most important 
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uncertainty being that in the low-energy region around the p° peak. The uncertainty in the 
hadronic vacuum polarization in this energy region may be reduced by combining the e+e~ 
annihilation data with those from p^v decay. There is also a hadronic contribution 

from light-by- light scattering diagrams, which has been estimated using chiral perturbation 



theory, and is thought to yield a smaller uncertainty in | |43|| . 

The estimate of the hadronic vacuum-polarization contributions [44| used in the E821 
paper does not include the latest e^e~ data from Novosibirsk |45] and Beijing [Q, nor 
the most recent r decay data from CLEO |^ 0. However, these are thought unlikely |^ 
to change the overall picture: we recall that the quoted hadronic error ~ 7 x 10^^'' is much 
smaller than the apparent discrepancy and the experimental error. Advocates of new physics 
beyond the Standard Model may therefore be encouraged. However, a final conclusion must 
await the publication of more — 2 data and the achievement of consensus on the hadronic 
contribution. 

A priori, the BNL measurement favours new physics at the TeV scale, and we consider 
the best motivated candidate to be supersymmetry. Even before the hierarchy motivation for 
supersymmetry emerged, the potential interest of was mentioned, and a pilot calculation 
performed ||9|, followed by many others 51]. Some time ago, it was emphasized |^ 



that the BNL experiment would be sensitive to a large range of the parameter space of the 
CMSSM with universal soft superymmetry-breaking parameters at the input GUT scale, 
determining in particular the sign of the Higgs mixing parameter /i |^. A large number 
of theoretical papers have discussed the interpretation of the BNL measurement within 
supersymmetry |5^, These calculations generally agree that /i > is favoured by the 
BNL measurement. The calculations we use in this paper are taken from which are 
based on |^ 0, including also the leading two-loop electroweak correction factor 



In this paper, we do not impose the BNL — 2 constraint in the form (|^), though we 
do bear it in mind in the selection of points, for example in the relative weighting of points 
with /i > and /i < 0. About half of the points we propose yield values of that are 
compatible with (0) within two standard deviations, and several of the points lie within one 
standard deviation. Table shows the values of 6a^ calculated in our proposed benchmark 
scenarios. 



^More data on r decays can be expected from the LEP experiments and the B factories. 
^For other recent calculations, see 1531 ■ 
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2.4 Cosmological Relic Density 

Like most analyses of CMSSM phenomenology for future colliders, we assume that R parity 
is conserved. This implies that the lightest supersymmetric particle (LSP) is stable, and 
hence should be present in the Universe today as a cosmological relic from the Big Bang, 
constituting part of the dark matter. If the LSP had either strong or electromagnetic inter- 
actions, it would bind with conventional matter to form anomalous heavy isotopes. These 
are not seen down to levels far below the calculated relic density, so the LSP can have only 
weak and gravitational interactions 0. There are scenarios in which the LSP is not the 
supersymmetric partner of any of the Standard Model particles. For example, it might be 
the gravitino or axino. In these cases, cosmological constraints on the dark matter density 
cannot be used to constrain the CMSSM in a useful way, and values of mi/2 and mo larger 
than those we discuss would also be allowed. 

Among the supersymmetric partners of Standard Model particles, LEP data and direct 
searches for the scattering of cold dark matter particles appear to exclude the possibility 
that the LSP is a sneutrino u in the MSSM The most viable LSP candidate seems to 
be the lightest neutralino aiid this is the hypothesis adopted here. Since the sparticle 
spectrum is explicitly calculable in the CMSSM, we concentrate on regions of its parameter 
space in which the LSP is a neutralino, to the exclusion of other regions. 

Astrophysics and cosmology provide many independent arguments that most of the grav- 
itating matter in the Universe is invisible. Some of this is certainly baryonic, but the con- 
sistency of cosmological nucleosynthesis calculations with the observed light element abun- 
dances suggest that most of the dark matter is non-baryonic . This conclusion has been 
reinforced by recent estimates of the cosmological baryon density Qf, based on microwave 



background data, which suggest |^ 



n^h' ~ 0.02, (3) 

with an error of about 20 %, where h is the present Hubble expansion rate in units of 
100 km/s/Mpc. The Hubble Key Project ||60[ and other measurements indicate that h"^ ~ 0.5, 



again with an error of about 20 %. The estimate (|^) is much smaller than the corresponding 
estimate of the overall matter density pU[ : 



rimh^ ~ 0.14 ±0.04. (4) 

We conclude that most of the matter in the Universe is in the form of non-baryonic dark 
matter, and hypothesize in this paper that it consists mainly of the lightest neutralino x- 
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For the purpose of this paper, we assume 

0.1 < n^h^ < 0.3. (5) 

The upper hmit being a conservative upper bound based only on the lower limit to the age 
of the Universe of 12 Gyr. Larger values of Q^h"^ would require values of mi/2 and mo larger 
than those we discuss, in general. Smaller values of Q^h'^, corresponding to smaller values 
of mi/2 and mo, are certainly possible, since it is quite possible that some of the cold dark 
matter might not consist of LSPs. Axions and ultraheavy metastable relic particles are other 
candidates that might contribute. However, allowing smaller values of f^^^^ would open up 
only a very small extra area of the (mi/2, mo) plane, as we see shortly. 

We base our relic density calculations on a recent analysis ||2^ using SSARD that extends 



previous results |^ to larger tan/3 > 20. We note here two important effects on the 



calculation of ^'^^^ were discussed in [22|, which are due to improvements of previous 



calculations of x — ^ coannihilations and direct-channel XX annihilations through the heavier 



neutral MSSM Higgs bosons H and A |18, 19, 20, 21, 22]. Both of these effects extend the 



region of CMSSM parameter space consistent with cosmology out to values of mo and mi/2 
that were larger than those found at smaller values of tan/5 |1T5| , |6l| . As we discuss later, 
good overall consistency was found |^ between these relic density calculations, the LEP and 



other sparticle mass limits, the LEP Higgs 'signal' and measurements of 6 ^ s'j, and also the 
recent BNL measurement of — 2 if /i > 0. There is also a region of the (mi/2, '^o) plane at 
relatively large values of mo, close to the higgsino LSP area, termed the 'focus-point' region. 
This is consistent with 6 — 57 for any tan /?, and may also be consistent with — 2 if tan (3 



is large and mi/2 is relatively small, according to the BMPZ code ||33| although not according 
to SSARD for the input parameter values used in Fig. ^ (see the discussion in Sec. H). 
Table ^ shows the values of Q^^h'^ calculated in our proposed benchmark scenarios. 

2.5 Electroweak and Cosmological Fine- Tuning 

Here we discuss two distinct issues: the fine-tuning of CMSSM parameters that is required 
to obtain the electroweak scale, and the sensitivity of the cosmological relic density to input 
parameters. 

As mentioned in the Introduction, the TeV mass scale for supersymmetry is largely moti- 
vated by the gauge hierarchy problem: how to make the small electroweak scale m^ <^ mp ~ 
10^^ GeV 'natural', without the need to fine-tune parameters at each order in perturbation 
theory 0. This is possible if the supersymmetric partners of the Standard Model particles 
weigh ^ 1 TeV, but the amount of fine-tuning of supersymmetric parameters required to 
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obtain the electroweak scale increases rapidly for sparticle masses ^ 1 TeV. In an attempt to 
quantify this, it was proposed to consider the logarithmic sensitivities of the electroweak 
scale to the supersymmetric model parameters ai. 

ai dmz 



In the CMSSM with universal soft supersymmetry-breaking parameters, the fundamental 
parameters include the common scalar mass mo, the common gaugino mass mi/2, the 
common trilinear parameter Aq at the GUT scale, the supersymmetric Higgs mass parameter 
yU at the GUT scale and the supersymmetry-breaking Higgs mass parameter B at the GUT 
scale. These are the fundamental dimensionful parameters which are expected to be directly 
related to the physics responsible for breaking the electroweak symmetry and generating the 
correct size for the electroweak scale. In view of the sensitivity of the electroweak scale to 
the top (and possibly the bottom) Yukawa coupling Xt{Xb), some (but not all) of the Yukawa 
couplings at the GUT scale are sometimes included among the fundamental parameters in 
(^) []. In what follows, we quote the values of A for both cases: first, considering the fine 
tuning only with respect to the dimensionful CMSSM parameters, and then including also 
the sensitivity to At. 

An analogous measure of the amount of fine-tuning needed to obtain in the CMSSM a 
relic density Q^h"^ in the range preferred by cosmology has been proposed recently |p4 |: 

In this case, we hold mz fixed, and therefore, the set of input parameters {ai} becomes 
{niQ, mi/2, Aq, tan P, sgn{^)}. The relic density is also quite sensitive to the values of the 
Standard Model parameters rrit and mi,. 

Table ^ shows the values of A and A^ calculated in our proposed benchmark scenarios. 
The first (second) row for A shows the electroweak fine-tuning without (with) Xt included 
among the Oj (the A;, dependence of A is relatively mild). The first (second) row for A^ 
shows the cosmological fine-tuning without (with) rrit and rrib included among the a^. The 
dependence of A^ on is significant at high tan (3, particularly for point K, L and M. We 
see from Table |^ that, as a rule, the electroweak fine-tuning roughly scales with mi/2 and is 
independent of mg, if only the sensitivity to the dimensionful parameters is considered - this 
is in essence the focus-point phenomenon [^, |6^ . On the other hand, A^ behaves similarly 
to A, except in the 'focus-point' and rapid-annihilation 'funnel' regions, where it is found 



that there is a strong sensitivity of Qh to the input parameters [24 



''For an extensive discussion of the philosophy behind these choices, see 
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We emphasize that fine tuning should not be confused with instabihty: the fact that 
CMSSM model parameters might need to be adjusted carefully in some cases does not mean 
that the resulting points in parameter space are inherently unstable. They are perfectly good 
electroweak vacua, and cannot be excluded a priori. The extent to which one cares about 
the amount of fine tuning depends on the underlying measure in CMSSM parameter space, 
which is of course unknown at present. Moreover, there are surely correlations between the 
input parameters, and it is known that these may reduce radically the apparent amounts 
of fine tuning. Finally, extending the spirit of @ and (^, one might define alternative (or 
additional) measures of fine-tuning refiecting the sensitivities of other physical observables 
to the fundamental parameters, which might change again the relative weights of the points. 



For example, one could consider large CP- violating phases |65] and ask about the degree of 
fine-tuning required to bring various electric dipole moments in accord with experiment. One 
would then find that in this sense the bulk points are much more fine-tuned than the focus 



points [^|, coannihilation 'tail' points |15] and rapid-annihilation 'funnel' points [^|. For 



these reasons, measures of fine tuning come with impressive health warnings on the packet. 
Hence, we do not use their values as selection criteria for benchmarks. However, we do see 
clearly that some models are more finely tuned than others. The values of A vary by a factor 
of 28 between benchmarks B and M (27 if is included among the Oj.), whereas the values 
of vary over a factor 81 between benchmarks B and E (840 between benchmarks B and 
F if rrit is included among the a^). 

2.6 Combination of Constraints 

The interplays of all these constraints in the (mi/2,'^o) planes for some values of tan (3 are 
illustrated in Fig. ^ and |^. The very dark (red) triangular regions at large mi/2 correspond 
to nif-^ < m^, where fi is the lighter f mass eigenstate. These regions are ruled out by the 
requirement that the LSP be neutral. We show as (red) dash-dotted lines the nih = 113 GeV 
contour calculated using FeynHiggs We see that the Higgs mass bound from LEP 

excludes regions of small mi/2, and has strongest impact at low tan/5. In a previous analysis, 
it was shown in that the whole of the plane which is of cosmological interest is excluded 



by the Higgs bound for values of tan/? ^ 3.5. Hence, for the benchmarks, only values of 
tan/? > 5 were considered. The (dashed) bound on the chargino mass from LEP excludes 
very low mi/2 values, almost independently of tan/5, and the LEP selectron constraint (dot- 
dashed) excludes a region around the origin in the {mi/2,mQ) plane. We do not show them 
on all the panels in Figs. |] and ||, but only in panels (a,b) of Fig.0 for tan/5 = 5, 10 and 
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/i > 0: their locations are similar for the other CMSSM cases studied. The branching ratio 
for b S'j excludes a dark (green) shaded area at low mi/2. Its impact increases with 
increasing tan/3, and is larger for /i < 0. 

The cosmological constraint 0.1 < ^^^"^ — ^-^ allows a region shown in light grey 
(turquoise), which exhibits a narrowing coannihilation strip that extends at large mi/2 into 
the domain where the fi is the LSP. This defines upper bounds on the allowed values of mi/2 
(and hence m^) in the coannihilation region, which are 

mi/2 ~ 1400 GeV, ~ 600 GeV, (8) 

for tan/? ^ 20, increasing at larger tan/?: mi/2 < 1900(2200) GeV is allowed for tan/? = 
35(50), as seen in (b, c, d) of Fig. |^. The 'tails' of these regions are potentially beyond the 
physics reach of the LHC, but are disfavoured by g^^ — 2. 

We also see in panels (c, d) of Fig. || the possibility of a 'funnel' extending to large mi/2 and 
mo where an acceptable relic density is made possible by rapid direct-channel annihilation 
XX ~^ H, A. There is also an allowed cosmological strip at large ttlq where the LSP has a 
significant higgsino component, and as a result, neutralino annihilation to gauge boson pairs, 
as well as s-channel Higgs exchange are enhanced [T7, 37|. This region lies in the 'focus-point' 



region and is present in all panels with tan (3 > 5, although for improved readability of the 
figures we choose to show it only in panel (d) in Fig. ^ Just above the focus-point region 
there is a shaded area with no acceptable electroweak symmetry-breaking solutions, and a 
light higgsino-like chargino near its boundary. Areas in Fig. ^ and § between the 'focus- 
point' and the other shaded cosmological regions have values of ^^^"^ that are too large, and 
hence are excluded by the cosmological relic density constraint. As already commented, the 
unshaded areas at lower mo values have Q^h"^ < 0.1, and hence are in principle allowed by 
cosmology^. However, the remaining parts of these regions compatible with other constraints 
are quite small |. 

Finally, we note that the — 2 result prefers the diagonal band at low mo and mi/2 
shown in darker grey (pink). The one-sigma band is indicated by dashed lines and the full 
lines represent the two-sigma band. We see that there is good overall compatibility between 
g^ — 2 and the other constraints for tan/? ^ 10 and /x > 0. The gfj, — 2 constraint disfavours 
large values of mi/2 and mo, excluding, for example, the tails of the cosmological region. 

^The central regions of the direct H, A annihilation channels in Figs. ||(c,d) are also allowed, as the relic 
density is very small there. However, the exact position of this region is sensitive to the input parameters, 
as the cosmological fine-tuning measure indicates. 

^We discuss this point again in Section 6. 
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3 Proposed Benchmark Points 



Supersymmetric benchmark points have a venerable history in physics studies for future 
coUiders |, |ll], H, |6| . They were useful in showing how many spectroscopic measurements 



might be possible at the LHC, and in demonstrating the precisions possible there and with 
an e~^e~ linear collider. However, only about 3 out of the 26 points previously studied 
are clearly compatible with all the LEP constraints, though some cases may survive if the 
theoretical errors in calculating rrih are favourable to them. The points used previously also 
could not take into account the recent constraints from b —>■ sj and — 2. Furthermore, 
many of the previous points also give unacceptably large relic densities. For instance, six 
points in CMSSM parameter space were studied in detail for the LHC None of them 
have survived the most recent LEP2 limits, dark matter constraints and Qfj. ~2 constraints. 

We have chosen our proposed new benchmark points for tan/? = 5, 10, 20, 35 and 50 to 
span the possibilities in the preferred regions. The locations of the points in the {1^1/2,^0) 
planes for different values of tan/3 are shown in Figs. |^ and 0. As already remarked, the 
points are probably all consistent with the LEP Higgs mass constraint mh > 113.5 GeV, once 
theoretical uncertainties are taken into account. We also took note of the (?^— 2 measurement, 
so that most points have /i > and several are within the 2-a experimental range. In some 
cases, points with different values of tan/3 give rise to very similar particle spectra and decay 
characteristics, and it was decided to keep only one example, so as to avoid duplication. 
For this and many other reasons, the chosen points should not be considered an unbiased 
statistical sampling of the CMSSM possibilities. However, we did make an effort to probe the 
different possibilities. Thus, we include two 'focus-point' models, two in the coannihilation 
tails at large mi/2, and two in rapid XX —>■ H,A annihilation funnels, and we kept two points 
with /i < 0. 

The input parameters of the benchmark points, labelled from A to M, are listed in Table |l[ 
We now make some comments on the individual points. 

A : The only allowed points for this small tan /3 = 5 are far into the coannihilation tail, 
and thus have relatively large mi/2, essentially to ensure rrih > 113 GeV. For this 
reason, this value of tan /3 is now disfavoured by (7^ — 2. It would be possible to choose 
a smaller value of mi/2 if one made greater allowance for theoretical error in the ruh 
calculation, e.g., by choosing nit > 175 GeV and Aq > 0, and relaxing the — 2 
constraint. We note that, when /i < 0, consistency with the Higgs limit would have 
required mi/2 > 830 GeV, and we do not consider this limited region for further study. 
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B : This point with tan/5 = 10 has much smaller mi/2, and hence mh < 113 GeV in our 
nominal FeynHiggs calculation. Though it formally fails the Higgs mass bound, it does 
so just barely: m/j = 112 GeV for this point, which is compatible with LEP [^] within 
the theoretical errors. On the other hand, it satisfies the g^~2 constraint within about 
one a. Thus points A and B take complementary points of view concerning these two 
constraints. 

C : This second point with tan/5 = 10 is compatible with — 2 at the 2-a level, as 
well as with rrih = 113 GeV in our nominal FeynHiggs calculation. It is therefore 
intermediate in philosophy between points A and B. 

D : Again with tan/? = 10, but one of just two points with < 0. This point has g^ — 2 
within about 3 cr of the current experimental value. Like point C, it is compatible with 
TTih = 113 GeV according to FeynHiggs. Note that the effect of the b ^ sj limit is 
similar to that of the Higgs limit over the range of mo favoured by the relic density. 
Both constraints force this point into the coannihilation region. 

E : The first of two 'focus-point' models with large mo at tan/5 = 10. Note that SSARD 
uses rrit = 171 GeV for this and the other focus point, whereas the ISASUGRA 7.51 
version uses m^ = 175 GeV. For this reason, FeynHiggs yields somewhat different 
values of m^ = 112, 116 GeV, respectively. Both of these are compatible with the LEP 
lower limit, taking into account theoretical uncertainties. With SSARD, focus-point 
solutions could also be obtained with the default choice of m^ = 175 GeV, but at 
larger mo for the same value of m,i/2- 

F : This second focus-point model has larger mi/2 and mo, and so has a higher Higgs 
mass: rrih = 115, 121 GeV for the SSARD and ISASUGRA 7.51 versions. Like point E, 
this point uses rrit = 171 GeV. Both points are about 2.5 a away from the central value 
of g^ — 2. We note that the focus-point region extends to larger mi/2 (and mo). 

G : One of two points with moderate tan /5 = 20, this one is chosen to have the relatively 
low value nih = 114 GeV. It is also just consistent with the b ^ s'y constraint, and 
agrees with the g^ — 2 measurement at the 1-a level. 

H : A second point with the same moderate tan/5 = 20, but this time at the end point 
of the coannihilation 'tail' with very large mi/2- Correspondingly, it has a larger 
Higgs mass: nih = 121 GeV, and a small value of g^ — 2, about 2.5 a from the central 
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experimental value. The small m^^ ~fT^x ^iciss difference has its own interesting features 
and challenges, as we discuss later. 

I : A very ((7^ — 2)-friendly point at large tan (3 = 35. Note that the inclusion of one-loop 
corrections to are essential here. If they are not included, this point has fi as the 
LSP. The 6 — » S7 constraint is the dominant one at small mi/2 for this value of tan/3. 

J : A second point with tan/3 = 35 and fi > 0, this time about half-way along the 
coannihilation 'tail' at large mi/2, that is compatible with (7^ — 2 at the 2-a level. 

K : One of two points in a rapid XX H,A annihilation 'funnel', and one of just two 
points with fi < 0. This point has tan/3 = 35, which is (almost) the largest value 
where we find consistent electroweak vacua for this sign of fi, with our default choices 
of the auxiliary parameters mt,mb and Aq. It is far from saturating the experimental 
constraints, apart from (7^ — 2. We note that the rapid-annihilation 'funnel' also extends 
to larger mi/2 (and itlq). 

L : One of two points with (almost) the largest value of tan /3 = 50 for which we find con- 
sistent electroweak vacua for /i > 0. This point lies in the ((7^ — 2)-friendly 'bulk' of the 
cosmological region, and is highly compatible with the other experimental constraints. 

M : A second point in a rapid xX ^ H,A annihilation 'funnel', again for (almost) the 
largest value of tan /3 = 50 allowed for fi > with our default choices of the auxiliary 
parameters. This point has small — 2, but satisfies all the other experimental 
constraints. 

It is characteristic of all the solutions that the lightest Higgs strongly resembles a Standard 
Model Higgs boson, whilst the other Higgses are heavy and nearly degenerate in mass. The 
LSP is in all cases almost a pure B, except in the focus-point cases, where a non-negligible 
H component is also present. 

4 Theoretical Uncertainties and Comparisons between 
Codes 

Our proposed benchmark points were chosen using the code SSARD to run all the input 
parameters down to the electroweak scale, impose the electroweak symmetry-breaking con- 
ditions there, and evaluate the physical chargino and neutralino masses including one-loop 
radiative corrections, which are generally C(5)%. This code does not include the one-loop 
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taking values that increase with 
33[]. As mentioned earher, the FeynHiggs code [BO] is used to 



corrections to slepton masses, which are generally 0{1 
tan j3 and decrease with 
evaluate the Higgs spectra. 

Some words of caution are in order. There are several available programs for calculating 
the supersymmetric particle spectrum and computing the physical observables considered 
in this paper. In general, these programs use different algorithms and, depending on their 
purpose, may have different levels of sophistication. The magnitudes of the resulting differ- 
ences depend on the particular observable considered, on whether the results are expressed as 
functions of GUT-scale input parameters or physical masses, etc.. For example, in Fig. ^ we 
show the cosmologically preferred region fl^ which was obtained with Neutdriver and 
mass spectra from the BMPZ code , for the same parameters as in Fig. We can see that 
the relic densities found for the same values of mi/2 and mg can be quite different, reflecting 
in part the sensitivity A*^ [24| noted earlier, which is generally larger at large tan/3 
In Figs. 0d and there are points at intermediate values of mo and mi/2 where the results 
for ^x^"^ '^^^ differ by as much as a factor of 20. We have sought to understand the origin 
of this apparent discrepancy and the related theoretical uncertainty in the calculations. 



m, = 175 GeV , tan p = 50 , |x > 
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m^ = 170.3 GeV , tan |3 = 52 , ^ > 
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Figure 4: Left: the cosmologically preferred region obtained with Neutdriver \7(\J and 
mass spectra from the BMPZ code j^^, for the same parameters as in Fig. ^d. Right: the 
corresponding result from SSARD, but for tan/3 = 52 and rrit = 170.3 GeV. 



Comparing codes, we find that the bulk of the effect is due to differences in the super- 



^°The results of [20| are not so dissimilar from those of |£2|. We are aware of other calculations underway 
(G. Belanger, F. Boudjema, A. Djouadi, M. Drees, A. Lahanas and L. Roszkowski, private communications), 
which also give rather diverse answers. 
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symmetric mass spectrum, most notably the values for the CP-odd Higgs mass and the 
fi parameter. We remind the reader that although both SSARD and BMPZ are NLO programs, 
numerical differences do arise at NNLO due to a different treatment of the NNLO terms. 
We list some of these effects below: 

• The treatment of gauge coupling unification. In SSARD, all three gauge couphngs are 
unified at the GUT scale, and the resulting prediction for the weak scale value of as 
is shown in Table |l|. On the other hand, in BMPZ (ISASUGRA), is treated as an 
input at the weak scale, with the following values: = 0.119 (0.118). The three 
gauge couplings do not in general unify, and model-dependent GUT-scale threshold 
corrections are assumed to account for the mismatch. This difference in the value of 
as affects the evolution of Xt between mz and mcuT, as well as the extraction of A;, 
from mf,{mi,). 

• Differences in the extraction of Xt{mz), which amount to a ~ 2.5% effect. Most of this 
uncertainty comes from using the running top quark mass (SSARD) versus the pole mass 
(BMPZ) as an argument in the one-loop correction to rrit 0. Given that the one-loop 
correction to is typically about 11%, the resulting differences in Aj are well within 
the NNLO uncertainty, but have a major impact (up to ~ 30% at large mo) on the 
extracted value of the n parameter. 

• Imposing electroweak symmetry breaking. One can choose to minimize the effec- 
tive potential at the scale mz (SSARD), or at the scale of the average stop mass 
Q = {TTif^^mf^}^^'^ (BMPZ and ISASUGRA). In the former case, care is taken to include 
corrections 0{ln{Q/mz)), in order to avoid spurious differences in the values of the 
/i parameter, which affects the gaugino-higgsino mixing and the pseudoscalar mass 
mA- The latter has a major effect in regions where s-channel annihilation through A 
exchange is dominant. 

• Computation of the Higgs sector. All of the effects already mentioned have an impact 
on the calculation of the Higgs boson masses. Furthermore, SSARD uses the results of 
p9|l , which include some two-loop effects, while BMPZ and ISASUGRA apply the one- 
loop Higgs mass corrections |Q at the scale Q. The tjia calculation is also affected by 
numerical differences in the extracted values of Xb- Altogether, we find that the values 
for rriA in the two cases can differ by as much as ~ 20%. 



"Both BMPZ and SSARD include the SM two- loop MS contribution 0, assuming it is close to the DR 
value. 
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Differences of such magnitudes are usually not important for collider phenomenology. 
However, this is not the case with the relic density calculation at large values of tan/3, 
because of the high sensitivity (^. Recall that then the dominant annihilation channel is 
s-channel A exchange into bottoms, and that the annihilation cross section scales as 

2 2 

^ ~ 2 fTf' 2,2 A^iA^3tan^/3, (9) 

where A^i^a are the Bino and Higgsino components of the LSP. The product NiN^ is relatively 
large when — Mif ^ M|. It is easy to see that in regions of parameter space where one 
is not very far away from the A pole, and yU is not very large, the ~ 20% effects in niA and 
/i can combine to an ~ 10 effect in the ^-xh"^ calculation. 

It is possible to vary slightly the input parameters and obtain reasonable agreement for 
the spectra obtained using SSARD and BMPZ/ISASUGRA 7.51. Panel (b) of Fig. ||b shows 
the cosmologically preferred dark matter region obtained with SSARD using slightly different 
inputs: tan/5 = 52 and = 170.3 GeV. We now see improved agreement (to within a 
factor of 2) between Figs, ^a and b. We have isolated several factors that contribute to this 
remaining difference, all of them related to the dark matter calculation. One is the treat- 
ment of the bottom-quark radiative corrections and bottom-quark mass, whose treatment 
in SSARD was discussed in |22[. Neutdriver typically uses a higher value of m^, which leads 
to a lower relic density. Neutdriver also does not integrate the Boltzmann equations, but 
uses an analytic approximation based on non-relativistic expansions of the annihilation 
cross sections. Other code differences concern the extent to which coannihilation effects are 
included, and the treatment of s-channel annihilation rates. Therefore, in comparing results 
over the CMSSM parameter plane, one must be sure not only to understand the differences 
in the spectrum codes, but also those in the codes used to calculate the relic density. 

The above discussion has focused on CMSSM spectra produced by two fully NLO codes, 
and we have seen that in spite of the present theoretical uncertainties, relatively well- 
matched spectra can be obtained, e.g., by varying the input parameters. Although very 
versatile, the codes generally available for Monte Carlo simulations of supersymmetry, such 
as ISASUGRA [||], SPYTHIA [|72[ and SUSYGEN [|73|, do not always include all of the ingredi- 
ents of a complete NLO analysis. For example, all three programs so far lack the one-loop 
radiative corrections to the chargino and neutralino masses 0. For the convenience of the 
experimental simulations that frequently use one of these codes, we have made searches in 
the input parameter space of ISASUGRA 7.51 to find points that reproduce specific features 

^^These are important not only for relating physical masses to GUT-scale input parameters, but also for 
relating searches for sparticle species, e.g., e+e^ X^X~ s-^id XiX2- 
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of the spectra of the different models shown in Table 1. In some cases, significant differences 
are inevitable, and compromises have been made. The non-implementation of radiative cor- 
rections to the chargino and neutralino masses is a particular snag, and we hope that this 
can be overcome in future issues of ISASUGRA, SPYTHIA and SUSYGEN. The ISASUGRA param- 
eters that best reproduce the most relevant features of the spectra of Table ^ are shown in 
Table ||, and we use them for the discussion of decay signatures in this Section. Despite these 
differences, we note that the general agreement between our calculated spectra and those 
generated by ISASUGRA 7.51 is goodj^. The most severe differences occur in the focus-point 
region and at large mi/2,'n^o at large tan/?. We note in passing some of the problems that 
arise when trying to match ISASUGRA 7.51 spectra. 

For this purpose, it is convenient to separate the proposed benchmarks into three classes 
of points. 

• Points in the 'bulk' and the coannihilation 'tail'. Here the most relevant masses to fit 
are those of the lightest neutralino x and the lightest stau fi, since the other particles 
are less relevant for the dark matter calculation. It is easy to fit and m^^ by varying 
just mi/2 and tjiq: one first varies mi/2 to fit the neutralino mass and then adjusts uiq 
to fit the stau mass. In these cases, the most relevant fact is that ISASUGRA 7.51 does 
not include the one-loop corrections to the neutralino and chargino masses, and hence 
returns lower neutralino masses for the same input values of mg and mi/2- In order to 
compensate for this, one has to crank mi/2 up in order that the tree- level mass from 
ISASUGRA 7.51 looks like the one-loop-corrected mass. However, the stop masses also 
depend strongly on mi/2, and therefore increase when this is done, as can be seen by 
comparing Tables ^ and |[ In turn, once the stop masses are higher, the Higgs mass 
also increases, as also seen in Table ^. It would, in principle, be possible also to vary 
the other input parameters so as to improve the match, but there is no well-defined 
procedure for doing this, and it is not clear what would be learned from such a lengthy 
exercise. Since the rest of the spectrum has intrinsic uncertainties anyway, we have 
not striven for perfect matches 0. 

• Focus points. These points are basically defined by the LSP mass and a certain value 

for the gaugino-higgsino mixing, which is in turn determined by the ratio of the one- 

""^^For example, for the same input parameters, the percent difference is generaUy not larger than 2% across 
the sparticle spectrmxi. The biggest differences (10-20%) are found for the heavy Higgs masses for points K, 
L and M. 

^''We note in passing that it was necessary to fix a minor glitch in ISASUGRA 7.51 in order to find a solution 
for point H. We are grateful to I. Hinchliffe and F. Paige for their help. 
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loop corrected values of the U{1) gaugino mass Mi, and fi. Here, the absence of 
one-loop corrections in ISASUGRA 7.51 is again unfortunate. In these cases, we have 
tried to match the masses of the LSP and X2> since they correspond to similar values 
of the effective Mi and fi values. Here again, we vary only a couple of parameters: 
mi/2 and mg, which determines fi and therefore the X2 mass. The change in mo alters 
the slepton masses by about 50 GeV, as seen by comparing Tables |l| and 

• Points in the rapid XX ^ A,H annihilation region. In these cases, the relevant masses 
are mi/2 and rriA- The first can still be fit, but rriA is only a very slowly varying function 
of mo. Since 171^ varies more rapidly with tan/5, we vary this parameter instead. For 
convenience, since we use a lower value of tan f3 for point M, we use it for point L as 
well. 

At the end of this lengthy discussion, it should be clear that the benchmark points are 
defined by the physical spectrum and not the values of mo, mi/2, tan/? etc., which are 
attached to them. The latter are nothing but convenient labels, which may vary from one 
program to another. In view of its versatility, in the rest of this paper, we use ISASUGRA 
to discuss sparticle decays and experimental signatures. The physics, of course, should 
remain unchanged, as long as the physical spectrum is the same. We stress that most of 
the differences discussed above are higher-order effects and represent in part the theoretical 
uncertainty in the calculation of the sparticle spectrum. 

5 Decay Branching Ratios 

The decay branching ratios in all the proposed scenarios have been evaluated using the pro- 
gram ISASUGRA 7.51. The benchmark points can be subdivided roughly into two classes, 
distinguished by the organization of their particle masses and their dominant decay modes. 
The most salient features of the decay signatures are displayed in Fig. |^, and can be sum- 
marized as follows. 

(i) For points A-D and G-M, we find that > M2, which has the following consequences: 

• The Xi* mass sequence starts as B^, with some admixture of Hf2- The remaining 
two states are mainly a combination of the two Higgsinos with some admixture of W^, 
a pattern typical of the gaugino region. Furthermore, mo < mi/2, giving at least one 
slepton lighter than the X2 xt^ and frequently several. 
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Supersymmetric spectra calculated using ISASUGRA 7.51 



Model 
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B 
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KQ) 


768 


343 


520 


662 


255 


548 


485 


1597 


454 


876 


1213 


560 


1842 




116 


113 


117 


117 


116 


121 


117 


124 


117 


121 


123 


118 


125 


RO 


893 


387 


584 


750 


1435 


2955 


521 


1813 


431 


851 


1070 


472 


1737 


AO 


891 


386 


583 


749 


1434 


2953 


521 


1812 


430 


851 


1069 


471 


1735 


H± 


895 


394 


589 


754 


1437 


2956 


527 


1815 


440 


856 


1074 


481 


1739 


X'l 


252 


98 


164 


221 


119 


434 


154 


664 


143 


321 


506 


188 


854 




467 


179 


303 


414 


197 


546 


285 


1217 


265 


594 


932 


349 


1558 


xt 


770 


349 


524 


667 


262 


551 


491 


1599 


460 


879 


1215 


564 


1843 


xl 
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540 


674 


317 


845 


506 


1608 


475 


889 


1225 


578 


1855 


xt 
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537 
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1217 
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xt 
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1225 
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1855 


9 


1357 
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932 
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804 


2372 


880 


3186 


828 


1669 


2516 


1051 


4029 




435 


206 


293 


383 


1433 


2942 


290 


1092 


308 


599 


1260 


450 


1957 


eR, fiR 


271 


145 


182 


239 


1427 


2897 


194 


709 


234 


425 


1088 


370 


1658 




428 


190 


282 


375 


1431 


2941 


278 


1089 


298 


593 


1258 


443 


1955 


n 


269 


137 


175 


233 


1415 


2873 


166 


664 


159 


334 


931 


242 


1249 




435 


209 


295 


384 


1427 


2930 


296 


1081 


319 


589 


1204 


439 


1809 




428 


189 


281 


374 


1425 


2929 


275 


1076 


285 


571 


1197 


409 


1803 


Ul, Cl 


1211 


546 


833 


1075 


1519 


3397 


789 


2834 


756 


1508 


2398 


978 


3789 


Ur, Cr 


1167 


529 


803 


1036 


1515 


3360 


764 


2716 


732 


1452 


2315 


948 


3643 


d-L, Sl 


1214 


552 


837 


1078 


1521 


3398 


793 


2835 


760 


1510 


2400 


982 


3790 


dR, Sr 


1161 


531 


801 


1032 


1515 


3356 


762 


2703 


730 


1445 


2305 


945 


3631 


tl 


940 


400 


635 


845 


987 


2401 


601 


2288 


569 


1190 


1883 


744 


3016 




1172 


580 


830 


1039 


1292 


2967 


785 


2649 


742 


1405 


2122 


918 


3378 


h 


1126 


503 


769 


998 


1281 


2961 


713 


2619 


647 


1335 


2053 


819 


3308 


62 


1161 


534 


803 


1028 


1503 


3333 


762 


2667 


725 


1406 


2121 


913 


3388 



Table 3: Mass spectra in GeV for CMSSM models calculated with ISASUGRA 7.51. The 
renormalization- group equations for the couplings and the soft superymmetry-breaking pa- 
rameters include two-loop effects, and the dominant one-loop supersymmetric threshold cor- 
rections to the third generation Yukawa couplings are included. The Higgs potential is mini- 
mized at the scale Q = {mi^mi^Y^'^. The Higgs and gluino masses are calculated at one loop. 
The rest of the superpartner spectrum is calculated at tree level at the scale Q. The input 
parameters have been adjusted so that the spectra best approximate those shown in Table [J. 
We have used the ISASUGRA 7.51 default values m^°^^ = h GeV and as{mz) = 0.118. It is 
assumed that Aq = 0, mj = 175 GeV. 
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Class (i) 



Class (II) 



Figure 5: Characteristic features of the spectra and principal decay modes in the two classes 
of benchmark points. 

• The X2 decays mostly to ii and uu and the xf to iu and ui. These decays are followed 
in most cases by £ — > £xi- 

• The X3 and Xi decay in 50 — 60% of the cases into xf + W^. Other principal decays 
are about 20 - 30% xl ^ X2^° and about 15 - 25% xl ^ X2^°- 

• The xf decays about equally into xt + xf + and xf + h^- 

• The squarks are in all cases heavier than any of the gauginos or sleptons and gluinos 
are heavier than squarks. They usually have large branching ratios for cascade decays. 

These benchmark points are therefore characterized by a large proportion of final states with 
leptons (plus jets). 

(ii) For points E and F, corresponding to the focus-point scenario, we find Mi < < M2, 
with the following consequences: 

• The neutralino mass sequence is roughly: x ~ X2 ~ X3 ~ Xa ~ but 
with large mixing among the states. As mo ^ rni/2, all sfermions are considerably 
heavier than the gauginos. 
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• The dominant neutralino decays are X2 3 Xi + Xs 4 + W"^, and to a smaller 

extent X4 ^ X2 + X3 + 

• For chargino decays, we find xf Xi + with essentially a 100% branching ratio, 
and xt ^ X2,3 + or xt + or xf + 

The dominance of decays involving or leads to final states with mainly jets and rarely 
leptons. These points all have a gluino lighter than the squarks. 

The above summary lists only the gross features of the benchmark points and, within the 
two classes of models, the individual points show important differences and will sometimes 
(e.g., solution K) deviate in significant details from the above statements. Taken together, 
the benchmarks cover a large variety of cases. 

Some examples of supersymmetric spectra in specific models are shown in Fig. Point 
C is in the 'bulk' of the cosmological region, in this case for tan/5 = 10 and /i > 0, point 
E is in the focus-point region at large mo, also for tan/? = 10 and /i > 0, point J is in 
the coannihilation 'tail' for tan/3 = 35 and yU > 0, and point M is in the rapid-annihilation 
'funnel' for tan/5 = 50 and /i > 0. Overviews of the dominant decay branching ratios of the 
various sparticles in these models are shown in Figs. |^ and ||. 

6 Prospective Supersymmetric Physics with Different 
Accelerators 

Estimates of the numbers of CMSSM particles accessible to different accelerators in the 
various proposed benchmark scenarios are summarized in Fig. ^. Caution should be used 
in the interpretation of this figure, because it does not capture the different qualities of the 
measurements at different colliders. For example, the masses and decay modes of weakly- 
interacting sparticles and Higgs bosons can be measured more precisely at e+e~ colliders 
than at the LHC, if they are kinematically accessible. Recall also that we have not chosen 
the points such that they give a fair representation of the more and less likely regions in the 
parameter space, but rather have opted to span the points over as large and diverse a region 
as possible, to allow studies of the consequences at the different colliders. Moreover, many 
of the sensitivities assumed require verification and refinement by future detailed studies. 
Hence Fig. ^ as such cannot be used to estimate if a particular collider would do well or 
not in discovering and measuring supersymmetry, but rather illustrates the diversity of a 
fraction of the possible scenarios and the complementarity of the different colliders. 
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Figure 6: The swper symmetric spectra and principal decay modes for benchmark points 
C,E,J and M, calculated usmg ISASUGRA 7.51. 
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Figure 7: Details of the principal decay branching ratios for sparticles in benchmark points 
C and E, calculated using ISASUGRA 7.51. 

The proposed benchmark points are ordered in Fig. ^ according to their degrees of com- 
patibihty with the recent measurement of — 2. Thus, if this measurement is confirmed 
as evidence for physics beyond the Standard Model, the reader can see immediately which 
benchmark scenarios are preferred, and how the prospects evolve for the different accelerators 
considered. One may also disfavour some of the benchmark points because of the amounts 
of hierarchical and/or cosmological fine-tuning they require, which we have documented in 
Table |. 

6.1 Tevatron 

The Fermilab Tevatron has just begun its next run, which is planned to deliver 1 fb^^ of 
data per experiment per year in its first two years, followed by a short shutdown for detector 
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Figure 8: Details of the principal decay branching ratios for sparticles in benchmark points 
J and M, calculated using ISASUGRA 7.51. 



maintenance and luminosity upgrade. In the subsequent years, the Tevatron experiments 
are hoping to collect as much as 5 fb~^ of data per experiment per year, which might be 
enough for a Higgs (or supersymmetry) discovery at the dawn of the LHC We now assess 
the prospects for these searches in the context of our proposed benchmark scenarios. 

• The search for the Higgs boson of the Standard Model is the cornerstone of the Tevatron 
Run II program. A considerable amount of effort has been put into optimising the Higgs 



discovery channels. The Higgs discovery reach in Run II is summarized in |£8[. Based 
on the current expectations for the performance of the new detectors, 3-cr evidence for 
(5-cr discovery of) a Standard Model Higgs boson with a mass of 125 GeV is possible 
with somewhat less that 10 (30) fb~^. Similar conclusions apply to the lightest CP- 
even Higgs boson of the MSSM, in the decoupling limit >> m^, and also for the 
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Figure 9: Summary of the prospective sensitivities of the various accelerators considered to 
CMSSM production in the proposed benchmark scenarios, which are ordered by their distance 
from the central value of — 2, as indicated by the pale (yellow) line in the second panel. 
We see clearly the complementarity between an e~^e~ collider and the LHC in the TeV range 
of energies, with the former excelling for non-strongly-interacting particles, and the LHC for 
strongly-interacting sparticles and their cascade decays. CLIC provides unparallelled physics 
reach for non- strongly-interacting sparticles, extending beyond the TeV scale. We recall that 
mass and coupling measurements at e~^e~ colliders are usually much cleaner and more precise 
than at hadron-hadron colliders such as the LHC. Note, in particular, that it is not known 
how to distinguish the light squark flavours at the LHC. 
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CMSSM. Therefore, according to the Tevatron study this machine will be able to 
discover the light Higgs boson in all of the benchmark points. 



In certain models with very light superpartners, the Tevatron also has a shot at finding 
supersymmetry [j73|. The gold-plated mode for supersymmetry discovery at the Teva- 



tron is the clean trilepton channel 3i ^t- The corresponding reach has been recently 
re-evaluated in |75|, with improved background estimates and optimized analysis 
cuts. In the case where the two-body decays of xf and X2 first generation sleptons 
are open, the 3-0" reach extends up to about mi/2 ~ 250 GeV. Point B therefore appears 
to be on the edge of the Tevatron sensitivity in this channel. However, by combining 
several additional channels, e.g., the like-sign di-lepton channel [^, the di-lepton 



plus tau jet channel |j78|, or the channels with jets, and isolated leptons 



as well as data from both collaborations, point B might be observable with the full 
data set. In all the other 12 cases, however, the superpartners are too heavy to be 
abundantly produced at the Tevatron, and will escape detection. We should note that 
to date there have been no dedicated studies of the Tevatron chargino/neutralino reach 
in the focus-point region, where both chargino states are often kinematically accessible. 

A very interesting case is illustrated by point H. The two lightest supersymmetric 
particles, fi and x?? are extremely degenerate Q, and the fi — Xi niass gap is smaller 
than the tau mass rrir. This possibility can occur anywhere along the borderline of 
the (red) shaded regions in Figs. || and ^, where m^^ — — 0^. In such a case, the 
two-body decay fi —>■ rxi is closed, and, in the absence of lepton-fiavor violation in the 
slepton sector, the light tau slepton predominantly decays via the four-body process 
fi iuiUrXi, and is stable on the scale of the size of the detector. In that case, by 
looking for long-lived massive charged particles 0, one can probe slepton masses 
up to 225 GeV. The specific example of point H falls beyond this projected sensitivity, 
but there exist points in the CMSSM parameter space where this signature is accessible 
to the Tevatron. 

The lighter stop, ti, is not within reach of the Tevatron collider for any of our proposed 
scenarios. However, we recall that we have always chosen Aq = 0, and that m^^ may 



be reduced significantly for some other values of Aq |83] 



We recall that the high level of degeneracy overcomes the kinematic suppression of the annihilation cross 
section due to the relatively heavy LSP mass. 
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This brief discussion shows that the Tevatron colhder has good prospects for discovering 
the hghtest CMSSM Higgs boson within all our proposed benchmark scenarios, and some 
prospects for detecting supersymmetric particles. In the latter searches, our benchmark 
points suggest some unconventional scenarios that should be kept in mind. 

6.2 LHC 

A preliminary inspection has been made of the LHC potential for these benchmark points, 
based on the simulation results summarized in the ATLAS Physics Technical Design Report 
and in the CMS Note A detailed study is clearly required before a real assessment 

of the LHC physics potential for these benchmarks can be made. For a preliminary look, 
the following assumptions were adopted to estimate the discovery potential of the LHC, 
assuming ATLAS+CMS combined, together with an integrated luminosity of 300 fb~^ per 
experiment. 

• The light Higgs is always within reach, since its production rate and decay signatures 
are very similar to those in the Standard Model, as demonstrated by detailed studies 
of Standard Model and MSSM light Higgs bosons at LHC. 

• The heavy neutral Higgses can be searched for via their decays H^, — > rr and /i/i, 
and the quoted if^ discovery range is based on the decay —>■ tv and ht. As is 
well known, there is a region at large itla and moderate tan (3 where the heavier Higgs 
bosons may escape detection at the LHC, and this 'hole' is reflected in our results 

• The observation of gauginos is either through direct Drell-Yan production of X2xf ? 
leading to trilepton final states, or via the inclusive production of charginos and neu- 
tralinos in the decays of squarks and gluinos. Neutralino decays to same-flavour di- 
leptons and missing energy yield a characteristic end point at the upper edge of the 
mass spectrum. If the leptonic decays of gauginos are enhanced because the sleptons 
are light, according to the published study |TU|, the 5 — cr discovery region for is 



bounded approximately by mi/2 < 200 GeV or mo < 0.45mi/2 for mi/2 < 900 GeV. 
However, it should be noticed that, if the squark and gluino masses are light enough to 
yield sufficient production rates, charginos and neutralinos with W- and Z-like branch- 
ing ratios are also observable from cascade decays (which applies to points E and L). 
The treatment of the major background to the clean trilepton channel in those studies 



-"^^Previous particle- level estimates of the LHC sensitivity can be traced from ||8J] . 
^^Supersymmetric decays of the heavier Higgs bosons might also be an interesting signature. 
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may now be updated along the lines of |76[. Heavier gaugino states are assumed 
observable in the same region of {1711/2,1710), provided their mass is less than 450 GeV. 
More detailed studies are needed to assess the observability of the various gauginos 
for the various benchmark points, since the results depend strongly on the involved 
masses and decay patterns. 

• The squarks have conservatively been considered observable when rriq < 2500 GeV 
and < 2000 GeV. However, we note that it is not known how to distinguish the 
different light squark flavours at the LHC, and that the separation of squarks from 
gluinos needs to be studied in more detail. 

• The charged sleptons can be observed through their direct production up to masses of 
~ 350 GeV and perhaps indirectly in the decays of other supersymmetric particles up 
to masses of ~ 250 GeV. In addition, we have assumed that they are observable only 
if the mass differences — m^o > 30 GeV. 

• As was mentioned earlier, the fi may be long-lived, if mi/2 and tuq have values close 
to the boundary where the fi becomes the LSP, as for point H. The LHC collabora- 
tions have already made analyses of their sensitivities for such scenarios, motivated by 
gauge-mediated models. They have found that such a fi can be detected at the LHC by 
measuring delayed signals produced by these particles in the external muon spectrom- 
eters of ATLAS 1^ and CMS whose time resolutions are ~ 1 ns. By combining 
time and momentum measurements, it is possible to determine nif^ . A mass resolution 
of ~ 3.5% was obtained for gauge-mediated point G2b in P], where m^^ ~ 100 GeV. 
At point H, the fi is much heavier, and its detectability requires further study. 

• The sneutrinos are not counted as observable, since they have large branching ratios 
for invisible decays into u uxi- 

Our preliminary estimates of the numbers of detectable particles of each kind for each 
benchmark point are summarized in Table ^. 

The great potential of the LHC for the discovery of squarks and the gluino as well as 
the lightest Higgs is clearly apparent. For most benchmark points, the discovery of some 
of the gauginos and sleptons would also be possible. However, it is also seen that for the 
points with the heaviest spectra, namely models F, H and M, the LHC may have difficulty 
in finding any MSSM state except the light Higgs boson. It is, however, not excluded that at 
least some of the states in these solutions will be found to be detectable after more detailed 
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Prospective observability at the LHC 
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Table 4: Numbers of particles for each benchmark model thought to be accessible at the LHC. 
The observabilities we assume are obtained by extrapolating from previous simulation studies 
by ATLAS and CMS. 



work has been done, as some squark masses are close to the limits chosen above. There is 
also the possibility of a luminosity upgrade for the LHC, which might extend its reach for 
squark and gluino masses up to about 3 TeV. We have not studied the extent to which the 
different squark flavours may be distinguished at the LHC. 



Previous studies [T^ demonstrated for six benchmark points that several sparticle spec- 
troscopic parameters could be measured with precisions of a few percent using kinematic 
distributions. This allows the fundamental parameters of the CMSSM to be constrained to 
1% — 10%, depending on the point studied. 

6.3 Linear e^e~ Colliders up to 1 TeV 

Electron-positron collisions at centre-of-mass energies up to about 1 TeV are envisaged by 



the TESLA NLC |8^, ^ and JLC 0] projects. In each case, a possible first phase 
at energies up to about 500 GeV is also considered. In both the lower- and higher-energy 
phases, data samples of the order of 1000 fb~^ or more will be collected over a period of 



several years. Recent supersymmetric benchmark studies for TESLA [jrT| have included two 
CMSSM scenarios with tan/5 = 3 (30), mo = 100(160) GeV, mi/2 = 200 GeV, fi > 0, and 
Ao= (600) GeV. In the case of the first point, the Higgs mass is less than 100 GeV, and 
thus now ruled out by LEP. 

Sparticles can be produced at any linear e+e^ collider if its centre-of-mass energy is 
larger than twice the mass of the sparticles, the pair production threshold, except for heav- 
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ier charginos and neutralinos, which can be produced in association with the hghtest chargino 
or neutrahnos, respectively 0. Typical supersymmetric signals are multi-lepton final states 
and multi-jet final states with large missing transverse energy. Sneutrinos can be detected 
at threshold energies if these can decay into channels including charged leptons with a suf- 
ficiently large branching ratio. For example, in some scenarios the Vj. can decay into tW . 
Otherwise, sneutrinos can be detected at higher energies via the decays of charginos. Apart 
from the threshold requirement to produce the particles, we require the branching ratio 
times cross section for the detectable channels to be larger than 0.1 fb in order to observe 
the sparticle, leading to at least 100 produced sparticles in the total data sample. 

Table |^ shows the sparticles which can be observed at a 1 TeV linear e+e~ collider 
in each of the different benchmark points proposed in this paper. The unpolarized cross 
sections and the decay branching ratios were computed using ISASUGRA 7.51. The listed 
numbers of observable particles take into account the decays of sneutrinos and neutralinos 
into undetectable final states. 

We note the following points concerning searches for CMSSM particles at a linear e^e~ 
collider in the energy range ^ 1 TeV. 

• The lightest CMSSM Higgs boson is always detectable in a first phase even below 
~ 500 GeV, and its mass, width, spin-parity and couplings can be measured with 



high precision |Tl|. In the second phase, interesting measurements could be made of 
the trilinear Higgs self-coupling. These possibilities would go far beyond the Higgs 
measurements possible at the Tevatron or the LHC, and could be used to constrain 
significantly the CMSSM parameter space. Higgs studies would form an attractive 
cornerstone of the first phase of a linear e^e~ collider. 

In all the proposed scenarios except H, K and M, a number of additional sparticles 
is within reach of a 1 TeV linear e^e~ collider. These exceptions are amongst those 
disfavoured by the measurement oi — In approximately half of the cases, partic- 
ularly those which are consistent with the g^ — 'i measurement at the 2 a level, a large 
fraction of the sleptons and gauginos can be detected. In particular, points B and I 
are very favourable for a linear collider with centre-of-mass energy up to 1 TeV. 

The LHC and such a linear collider are largely complementary in their sparticle mass 
reaches. For example, in the case of point L, the linear collider would see all the 
sleptons, while the LHC would see all the squarks. 



^®For sufRciently light neutralinos and sneutrinos, observation of the radiative production of otherwise 
invisible final states may be experimentally accessible. 
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Observable particles at linear e^e colliders 
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Table 5: Numbers of particles accessible for each benchmark model for various lepton- 
antilepton collider centre- of-mass energies in TeV. Channels are considered observable when 
their cross section times branching ratio to visible final states exceeds 0. 1 fb, taking account 
of the invisible final states originating from some neutralino and sneutrino decay modes. No 
considerations of realistic detection efficiencies have been included. 



In contrast to the LHC, one notices that the squarks are not generally accessible in 
the proposed models. 

Beyond the discovery of sparticles, a crucial issue in the understanding of the nature of 
any new physics observed will be the accuracy obtainable in the determination of the 
sparticle masses and decays, and also their quantum numbers and mixing. A strong 
advantage of lepton colliders is the precision with which such sparticle properties can 
be measured. Typically, the masses of sleptons and gauginos can be determined with 
a precision of a few per mille, by threshold scans and by measuring end points of two- 
body decay channel signatures in inclusive distributions. This uncanny precision, even 
for a limited number of sparticles, will be of cardinal importance for the reconstruction 



of the underlying supersymmetric model and breaking mechanism ||91|| . 
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Furthermore, the availabihty of polarized beams at a hnear colhder will provide addi- 
tional tools for identifying supersymmetric particles and allow for additional measure- 
ments of parameters of the supersymmetric model, such as the mixing angles of the 
sparticles ^2 . 



In a few cases the heavy Higgs particles are within reach of a 1 TeV e^e^ collider, 
but we also note that the reach for heavy Higgses can be extended with a photon 
collider option ||93|, in which the lepton beams are converted into high-energy photon 



beams, with a peak energy of up to 80% of the original incoming lepton beam energy. 
Heavy Higgses can be singly produced in two photon interactions, increasing the reach 
in detection to Higgs particles with masses of approximately 75-80% of the centre-of- 
mass energy. Such an option would be very useful for the points C, D, G and L where 
the heavy Higgs particles become accessible in the two-photon collider mode. 

• We note also that a precise measurement of the two-photon width of the light Higgs 
boson, which appears feasible at a photon collider |Q, might give evidence of the 
existence of new physics, such as super symmetry, even if no other new particles besides 
the Higgs boson have been discovered at the LHC or a 1 TeV linear collider. 

6.4 CLIC 

CLIC is a project for an e~^e~ linear collider which aims at a centre-of-mass energy of 3 
TeV, upgradable to 5 TeV at a later stage. The luminosity will be about 10^^cm~^s~^. To 



achieve this luminosity, CLIC will work in the high beamstrahlung regime ||95|. The physics 
potential of this machine is currently being studied |^ . 



To obtain a first, very crude, estimate of the potential of a linear e~^e~ collider in the 

multi-TeV range, we have used the same criteria as in the previous section, i.e., we assume 

that sparticles can be detected provided their production cross section multiplied by their 

branching ratio to visible final states exceeds 0.1 fb |^. CLIC is estimated to yield an 

integrated luminosity of about 1000 fb^^ in a year, of which about a third would be close 

to the nominal Ecm- Therefore, one could accumulate 1000 fb~^ at this nominal energy in 

about three years, corresponding to 100 events produced for a cross section of 0.1 fb. More 

accurate estimates of the CLIC sensitivity will become possible only after simulation of the 

signals and backgrounds. Preliminary studies indicate that, although CLIC will operate with 

^^This limiting cross section is why some sparticles are not counted as observable, even though they are 
kinematically accessible. For example, at 5 TeV , at point K the d/j and sr have cross sections of only 0.01 
fb and the 62 only 0.04 fb, whilst at point H the ii is produced with a cross section of only 0.03 fb. 
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more beamstrahlung and hence a less well-defined centre-of-mass energy than lower-energy 
linear colliders, there are several physics topics for which the accelerator environment is not 



a serious disadvantage 96 



The unpolarized cross sections and the decay branching ratios were again computed using 
ISASUGRA 7.51. The numbers of observable particles, after taking into account the decays 
of sneutrinos and neutralinos into undetectable final states, are also summarized in Table |[ 
In preparing this Table, the possibility of detecting sneutrinos from the two-body decays of 
charginos has been taken into account, and we note that, in scenario B, the sneutrinos are 
only observable through their associated production. Note that due to the decrease of the 
cross section as function of energy, the and i)r which are observable at 3 TeV but are 
in principle no longer at an energy at 5 TeV if the luminosity is the same as for a 3 TeV 
collider. However this would ignore the fact that the machine could run at lower energy, 
and also the fact that CLIC makes an 'autoscan', by virtue of the beamstrahlung photon 
spectrum. Hence, for the counting in Table |^, these particles are considered to be detectable 
at a collider which can go up to 5 TeV. Finally, we note that, at point E, the gluinos are 
observable at energies > 3 TeV, as they constitute the dominant decay mode of the squarks, 
but they are not listed in the Table. 

Assuming that the LHC and an e^e~ linear collider in the range ^ 1 TeV have been 
taking data for several years before the start of a 3 TeV machine like CLIC (CLIC3000), we 
infer that supersymmetry will most probably already have been discovered by that time, if it 
exists. Hence the role of CLIC may consist mainly of completing the sparticle spectrum, and 
disentangling and measuring more precisely the properties of sparticles already observed at 
the LHC and/or a lower-energy e^e~ linear collider. However, a machine like CLIC would 
be needed even for the direct discovery of supersymmetry in the most problematic cases, 
namely scenarios H and M. 

A few benchmark points emerge as typical of situations which could arise in the future. 

• Point C has very low masses, and is representative also of points A, B, D, G, I, L. 
In these cases, LHC would have discovered the if^, as well as seen the /i^, and also 
the gauginos x?, X2 ^^"^ Xi-, some of the charged sleptons, the squarks and the gluino. 
A 1 TeV linear collider would enable the detailed study of the and of the same 
gauginos and sleptons, and it might discover the missing sleptons and gauginos in 
some of the scenarios. However, one would require CLIC, perhaps running around 2 
TeV, to complete the particle spectrum by discovering and studying the heavy Higgses 
and the missing gauginos. CLIC could also measure more precisely the squarks and 



40 



in particular disentangle the left- and right-handed states and, to some extent, the 
different light squark flavours. 

Point J features intermediate masses, rather similar to point K. Here, the LHC would 
have discovered all the Higgs bosons, the squarks and the gluino, but no gauginos 
nor slcptons. The 1 TeV e+e^ linear collider would study in detail the and could 
discover the e^j, fi^ and fi, but other sparticles would remain beyond its kinematic 
reach. CL1C3000 could then study in detail the heavy Higgscs. It would also discover 
and study the gauginos and the missing slcptons, and even observe in more detail a 
few of the lighter squarks that had already been discovered at LHC. However, to see 
the remaining squarks at a linear collider would require CLIC to reach slightly more 
than 3 TeV. 

Point E has quite distinctive decay characteristics, due to the existence of heavy slep- 
tons and squarks. In this situation, the LHC would have discovered the all squarks 
and the gluino. The gauginos are in principle accessible, but their discovery may be 
made more difficult because their predominant decays are into jets, contrary to the 
previous benchmark points, and sleptons would remain unobserved. At a 1 TeV e'^e~ 
linear coUider, the detailed study of the hP and of the gauginos could be undertaken. 
The discovery of the flrst slepton, actually a z/g, could be made at CLIC3000, which 
could also study the three hghtest squarks. The discovery and analysis of the heavy 
Higgses would then require the CLIC energy to reach about 3.5 TeV, which would also 
allow the discovery of all sleptons and the observation of all squarks. 

Point M has quite heavy states, rather close to scenario H. The LHC would only discover 
the h^, all other states being beyond its reach. The existence of supersymmetry might 
remain an open question! The 1 TeV e+e~ linear collider also sees only the h^, but can 
study it in detail, and so might provide indirect hints of supersymmetry. A CL1C3000 
would be able to discover most of the gauginos and the fi, but CLIC with 3.5 to 4 TeV 
would be required to discover the remaining gauginos and sleptons. To discover the 
squarks, £^£^ collisions above 6.5 TeV would be needed. There is currently no e+e~ 
project aiming at such energies, and neutrino radiation becomes a hazard for lJi^lJi~ 
colliders at such energies. 

Point F also has heavy states. Here again, LHC sees only the hP, but would also find 
the gluino. The hghtest gauginos would be within the reach of a 1 TeV e^e" linear 
collider, and the heavier ones within the reach of CLIC3000. But it would only be with 
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collisions in the 6.5-7.5 TeV region that the heavy Higgses, the sleptons and the 
squarks could be found. 

• As in the 1 TeV e^e~ linear collider photon collider option for CLIC would 
extend the discovery range for heavy Higgs bosons. It would additionally allow one to 
discover all four Higgs bosons in scenarios E, H and M, for a 3 TeV collider, and also 
in F, for a 5 TeV Collider. 

• We also note that polarization would have advantages similar to those discussed earlier 
for a lower-energy e^e~ linear collider. 



7 Conclusions and Prospects 

We have proposed some benchmarks that span the possibilities still allowed in the CMSSM, 
following the explorations made by LEP. A grand summary of the reaches of the various 
accelerators is presented graphically in Fig. The different levels of shading (color) present 
the different types of sparticle: Higgses, charginos and neutralinos, sleptons, squarks and 
gluino. The first six points (I, L, B, G, C, J) are presently favoured: they are compatible 
within 2 a with the present — 2 measurement, and the fine tuning is relatively small for 
most of these points. Fig. ^ summarizes the discussion of this paper and exposes clearly 
the complementarity of hadron and electron machines. It is apparent that many alternative 
scenarios need to be kept in mind. 

Beyond the CMSSM framework we have discussed in this paper, one should consider more 
general versions of the MSSM, in which the GUT universality assumption is relaxed [0. We 



note that the phenomenology of gauge-mediated models [|T2| with long-lived neutralino NLSP 
is very well covered by our suggested CMSSM benchmark points, whilst a charged long-lived 
NLSP is found for our point H. Gauge-mediated models with a promptly-decaying NLSP give 
rise to distinct signatures, but there are very few theoretical models that predict these. The 



minimal gaugino-mediated models |jT3[ predict spectra which are very similar to those 



of the benchmark points in the 'bulk' region discussed above, and are therefore to a large 



extent covered by our analysis. Anomaly mediation provides an interesting framework |14 



for model building A 'phenomenological' model of anomaly mediation \\10(\\ has been 



incorporated in ISASUGRA, and might serve as a basis for future studies. Finally, we recall 
that R-violating models generally do not provide a suitable dark matter candidate, and 
contain so many new /^-violating couplings that all the possibilities cannot be covered with 
only a few benchmark choices. 
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As we have discussed, the Tevatron has a chance to make the first inroads into the 
spectroscopy of the CMSSM models we have studied, with the best chance being offered by 
the hghtest Higgs boson, followed by chargino/neutralino searches. 

The LHC is expected to observe at least one CMSSM Higgs boson in all possible scenarios, 
and will in addition discover supersymmetry in most of the models studied. However, we do 
observe that the discovery of supersymmetry at the LHC is apparently not guaranteed, as 
exemplified by benchmarks H and M. It would be valuable to explore the extent to which 
precision measurements at the LHC could find indirect evidence for supersymmetry in such 
scenarios. We note also that, in these cases, the squarks and gluinos he not far beyond the 
nominal physics reach of the LHC, and an upgrade of the luminosity to 10^^ cm~^s~^ might 
bring them within reach. 

An e"*"e~ linear collider in the TeV range would in most cases bring important additional 
discoveries, exceptions being benchmarks H and M, and possibly E. Moreover, such a linear 
collider would also provide many high-precision measurements of the Higgs boson and super- 
symmetric particle masses and decay modes, that would play a pivotal role in first checking 
the CMSSM assumptions and subsequently pinning down its parameters, or those of the 
model that supplants it. As such, it will be an essential tool for securing the supersymmetric 
revolution we anticipate. 

In many of the scenarios proposed, the discovery and exploration of the complete set of 
supersymmetric particles, and especially some of the heavy Higgses, gauginos and sleptons, 
will have to await the advent of a machine like CLIC, which may need to run at an energy 
considerably higher than 3 TeV. In particular, points F and M are very challenging and 
would need a coUider with a centre-of-mass energy of up to 7.5 TeV in order to discover 
all the CMSSM particles. Distinguishing the different squark flavours could be an interesting 
challenge for CLIC. 

If the CLIC technology cannot be extended this far, here might come a role for a high- 
energy ^1'^ ^1" collider, if its neutrino radiation problems could be overcome. A lower-energy 
collider would also be interesting, with its unique capabilities for Higgs measurements. 

History reminds us that benchmarks have a limited shelf-life: at most one of them can be 
correct, and most probably none. In the near future, the CMSSM parameter space will be 
coming under increasing pressure from improved measurements of g^^u — 2, assuming that the 
present theoretical understanding can also be improved, and h — > 57, where the B factories 
will soon be dominating the measurements. We also anticipate significant improvement in 
the sensitivity of searches for supersymmetric dark matter. 

We also note that astrophysical and cosmological estimates of the cold dark matter 
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density are converging, which may improve the upper hmit on Vt^Ji^. As we mentioned earher, 
lower values of Q^h!^, even below 0.1, would be possible if there exists some additional form 
of cold dark matter. In this case, the sparticle spectrum might be somewhat lighter than 
in the benchmarks we propose, though the scope for this is rather limited. An interesting 
question for the future is whether the sparticle measurements at the LHC and particularly 



a linear e~^e collider will enable accurate calculations of Q^h"^ to be made | 101 , 24]. This 



certainly seems possible in many of the benchmarks we propose, though probably not for 
those with very large values of A^. 

Needless to say, the preliminary observations presented above need to be confirmed by 
more detailed exploration of the benchmark scenarios we propose, and of any others proposed 
within the context of different supersymmetric model assumptions. These more detailed 
studies would certainly benefit from improvements in the available simulation codes. All of 
the necessary ingredients for a complete NLO analysis are already available |^|. In order to 
achieve the necessary precision, one now has to include the complete one-loop corrections to 
the physical particle masses as well as perform a careful treatment of sparticle thresholds, 
the effective potential and the Higgs boson mass parameters. 
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